clc,clear,close all
load('pic1_1_cz.mat')
%由于原图已经黑白分明所以不需要亮度变换和滤波
i=imbinarize(z);
% se = strel('square',3);
% se = strel([0 1 0;1 1 1;0 1 0]);
se = strel('square',30);
imshow(i)
for n=1:5
i= imdilate(i,se);
end
 for n=1:5
i= imerode(i,se);
 end
 figure
 imshow(i)
i=edge(i,'prewitt');%使用prewitt算子进行边缘检测
figure;
% imshow(i);